Order and mobility of solid vortex matter in oscillatory driving currents 
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We study numerically the evolution of the degree order and mobility of the vortex lattice under 
steady and oscillating applied forces. We show that the oscillatory motion of vortices can favor an 
ordered structure, even when the motion of the vortices is plastic when the same force is applied in 
a constant way. Our results relate the spatial order of the vortex lattice with its mobility and they 
are in agreement with recent experiments. We predict that, in oscillating applied forces, the lattice 
orients with a principal axis perpendicular to the direction of motion. 
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Two limiting cases have been distinguished in the de- 
pinning process of a steady driven vortex lattice (VL) 
in type-II superconductors. For a high density of weak 
short-range pinning centers, the motion of the VL is in- 
homogeneous only in a narrow region near the critical 
force, F c (which is determined by the collective pinning 
theory Q), whereas for strong pinning centers (or a small 
concentration of them) plastic deformations become im- 
portant and the motion is disordered j| g]. In the last 
case, the VL orders and shows elastic flow only for high 
enough driving forces ^, ||, |) (above a threshold force, 
F T ). 

Recently, there has been much interest in the evolu- 
tion of the order and mobility of the VL for different 
thermomagnetic histories and different temporal depen- 
dencies of the driving force, F L , which is provided by 
a flowing current. This is motivated by the observation 
of a variety of puzzling phenomena that include history- 
dependent dynamic response and memory effects [Q, |^, [)j . 
The nature of the ac response observed in transport and 
ac susceptibility measurements has made clear that it is 
not possible to directly extrapolate the force-dependent 
evolution of a steady driven VL to explain the new exper- 
imental findings for oscillatory forces. Some results are 
intrinsic to the oscillatory dynamics. In particular, re- 
cent experiments suggest that while for a given intensity 
of the dc current the VL disorders, an ac current of the 
same magnitude assists the VL in ordering J?], g|. This 
conclusion, however, is indirectly drawn from changes 
in the VL mobility. A key issue that naturally surges 
is the microscopic dynamics of vortices and the relation 
between VL's order and mobility with different temporal 
evolutions of the driving force. Nevertheless, in spite of 
the growing amount of experimental results on the ac dy- 
namics of vortices, this relation has not been investigated 
numerically up to now. 

In this Letter, we present a numerical study that ad- 
dresses simultaneously the order and mobility of the VL 
in ac and dc driving forces. We show that many of the 
experimental observations can be reproduced for a VL 



around the crossover between plastic and elastic regimes. 
For F L < F c , the VL is at rest whereas for F L > F T 
the vortex motion is ordered. For F c < F L < Ft, if 
the force is applied steadily, the vortex motion is always 
plastic and disordered but, if it is applied in an alternat- 
ing way, it can bring the VL to a more ordered state. 
We also find that a lower density of defects in the VL 
is concomitant with a higher mobility and that the VL 
tends to order with its principal lattice vector aligned 
perpendicular to the direction of motion in contrast with 
the orderly array of vortices in the steady force case . 
All of these observations may have important implica- 
tions on the interpretation of recent neutron scattering 
experiments 
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We consider a transverse 2D slice (in the x — y plane) 
of an infinite superconducting slab that contains rigid 
vortices that are parallel to the sample surface (H=_Hz) . 
The overdamped equation of motion of a vortex in posi- 



tionr, isF, = E™^ F w (r<->"i)+Efc* F^(r l -r£)+F i + 
Fj = V^t wnere F« is the total force on vortex i due to 
vortex- vortex interactions (F vv ), pinning centers (F vp ), 
the driving current J (F L ~ </>o J x z) and thermal fluctua- 
tions (Ff ). Here, ij is the Bardeen-Stephen friction coef- 
ficient, N v the number of vortices, N p the number of pin- 
ning sites and the location of the fcth pinning center. 
The vortex-vortex interaction is modelled by a modified 



Bessel function F™ (r^ — rj 



to . 



A V K 1 (|ri-rj|/A)fij 



where the dimcnsionless prefactor, A v , tunes the rigidity 
of the VL (|]. We cut off this interaction at 4A. Quenched 
random disorder is modelled by pinning centers in uncor- 
related random positions that exert an attractive force 



on vortices: F vp (ri — r?) = — A p 



T^. A p tunes 

the strength of this force and r p its range. The ther- 
mal force was implemented with a Box-Muller random 
number generator and has the properties (F?) — and 
(FT(t)FT(t')) = 2r)k B T5 ij 5 fi J(t - t') for a given tem- 
perature T. We normalize length scales by A and forces 

j 2 

by /o = ■ We consider a sample with periodic 

boundary conditions in the x — y plane and size Lx x Ly 
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with L y /L x — v3/2. The normalized vortex density is 
n v = N V X 2 1 L x L y — B\ 2 /(f)Q. We choose n v = 1 and 
a much higher density of pinning centers n p = 25 with 
r p = 0.2A so that the resulting random pinning poten- 
tial varies smoothly on a length scale of the order of, or 
higher than, the coherence length, £ (< A/10). A p has 
a Gaussian distribution with central value A r p = 0.02/o 
and a standard deviation of O.IA™. System sizes from 
400 to 1600 vortices have been investigated. We show 
results for 1600 vortices and 4 x 10 4 pinning sites and we 
set T ~ (no thermal noise) unless noted to the contrary. 

History effects have been frequently associated with 
changes in the density of defects in the VL. If this were 
the case, information about the history of the sample 
could not be retained either by a perfectly ordered VL 
or by a completely amorphous one, i.e. it would be re- 
quired a balance between the elastic and pinning ener- 
gies. In dynamical terms, this implies that, in samples 
that present a path dependent response, the VL would 
be in the region of the plastic to elastic crossover |2|, |). 
We start by identifying this crossover by looking at the 
shape of the velocity- force dependence (V — I), the vor- 
tex trajectories and their translational order for differ- 
ent degrees of softness of the VL which are attained by 
changing A v ||. Moving vortices induce a total electric 
field E = f v x z, with v = -|- £\ $j±. To quantify the 
degree of order of the VL we compute the structure fac- 
tor, S(k) = e jk ri |, and determine the average 
concentration of vortices with coordination number not 
equal to 6, n^efi using the Delaunay triangulation proce- 
dure. As stated in Ref. ||, in the plastic regime there is 
a marked upward curvature in the V — I curve which is 
not present in the elastic region. In addition, the thresh- 
old force Ft diminishes with A v and when approaching 
to the elastic regime it mingles with F c . This provides 
a simple criterion to identify the crossover that, in our 
case, takes place at A v ~ 0.85 — 0.9. Following this result, 
we perform the simulations at fixed A v = 0.8, which is in 
the region of the crossover but in the plastic side || [l2| . 

In order to study the oscillatory dynamics and, at 
the same time, to be able to directly compare the ob- 
tained results with those of steady forces, we choose 
square oscillating forces of strength F L and frequency 
u> = 2tt/P (F l = F L x). For steady driving forces be- 
tween F c ~ 0.0085/o and F T = 0.012/ the vortex flow is 
always disordered with a high density of defects, ridef (of 
the order of 20%). For higher driving forces, ridef dimin- 
ishes and the movement is ordered with all of the vortices 
moving at the same average velocity || . Experimental 
results U suggest that when vortices perform an oscilla- 
tory motion whose amplitude is of the order of the lattice 
constant, the healing of defects should be important. For 
this reason, we initially choose F L = 0.0118/o and P so 
that F L P/4 ~ 1 (in normalized units). 

Some of the main results of our investigation are pre- 
sented in Fig. [I] and ^. Fig. |l| clearly pictures the healing 



of defects that results from the oscillatory motion of vor- 
tices. We start with a perfectly ordered VL. Then we 
apply a steady force F L = O.OII8/0 m the horizontal 
direction. The left panel shows a snapshot of the con- 
figuration of the VL after reaching a stationary state in 
which the average velocity of vortices and the average 
concentration of defects remained constant. This state 
presents a relatively high concentration of defects. We 
then turn to an alternating driving force which reverses 
periodically but keeps the same strength (O.OII8/0). This 
alternating force is applied to the disordered state in the 
left panel. After 500 cycles, it obtained the configuration 
of the VL shown in the right panel. The corresponding 
structure factors of both configurations are also shown. 
The Delaunay triangulation and the structure factor re- 
veal an important reduction in ridef when the oscillatory 
force is applied. It is fundamental to emphasize that the 
initial state was obtained with a steady force with the 
same strength starting with a perfectly ordered VL. This 
implies that, for certain values of the applied force, the 
flow of vortices can be plastic and disordered introducing 
defects in the VL but, if the same force is applied in an 
alternating way, the defects heal and the VL reorders. 

In the Fig. ||a), we show the evolution of this reorder- 
ing. We plot the concentration of defects, ridef > the aver- 
age value of the absolute velocity, (V), and the average 
quadratic displacements of vortices per cycle, (A 2 xn) 
and (A 2 ypf), as a function of the number of cycles, 
N. We define (A 2 ^) = ^Ef"M^) " ^{N - l)] 2 
and (AV) = -fEf v MN) ~ UN - l)] 2 , where 
x t {N - 1) = Xi(N - 1) + Al cm (iV) and y,{N - 1) = 
yi (N-l)+AY cm (N) with (AX cm (N), AY cm (N)) a cor- 
rection due to the displacement of the center of mass in 
the cycle N. The results shown correspond to the average 
over 4 different random distributions of the pinning cen- 
ters. We see that, for an increasing number of cycles, ridef 
clearly decreases and the vortex mobility is enhanced (as 
reflected by a 40% increase in the average vortex veloc- 
ity). We also note that (A 2 xn, Un) diminish and that all 
magnitudes vary in an approximately logarithmic way as 
observed experimentally 0, |j| . The combination of these 
observations indicates that, as the number of defects de- 
creases, vortices are clearly more mobile and that, af- 
ter completing one oscillation, they return closer to their 
original positions ((A 2 xn) and (A 2 z/at) decrease). This 
suggests that vortices organize and move more coherently 
as they are forced to oscillate. 

In Fig. |^b), we show what happens when we vary 
P(F L ) while keeping F L (P) fixed. We plot ridef as a 
function of the parameter F L P/4 after applying 100 cy- 
cles of the square force to an initial disordered state as 
in Fig. ^] and Fig. ||a). An important result is that, if 
the period P is increased so that the excursion of vortices 
greatly exceeds the lattice constant, the VL distorts and 
does not reorder. This is not surprising because if the 
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FIG. 1: Delauney triangulation and structure factor, S(k) 
for steady forces (left) and oscillatory square-forces (right). 
The more-ordered vortex configuration at the right panel was 
obtained after 500 cycles of the oscillatory force applied to 
the disordered state in the left panel. See text. A v = 0.8 and 
Fl = 0.0118/q (in the horizontal direction). 



amplitude of the oscillation is high enough, the system 
should behave in the same way as when driven by steady 
forces. This is a prediction that should be observed in 
experiments as those in Ref. ^[ For low enough P, vor- 
tices oscillate in their pinning sites and the VL does not 
reorder either. 

Notably, if there is a tiny asymmetry in the ampli- 
tude of the square force (~ 5%) the VL quickly disorders 
and, after a few cycles, both ridef and the vortex mobility 
reach values close to those found with steady forces. This 
observation is reminiscent to the experimental results re- 
ported in Ref. |^ where the vortex lattice is observed to 
quickly become less-mobile when shaken by a sawtooth 
magnetic field. The reduction in the mobility of the VL 
has been attributed to disorder produced by a ratchet- 
like tearing of the VL due to the different Lorentz forces 
involved in the ramp up and down of the field ||. We 
should note, however, that in our simulations there is a 
net displacement of vortices after completing one cycle of 
the asymmetric waveform that is not expected to occur 
in the experiments. 

An interesting issue surges in relation to the orienta- 
tion of the VL. At high-enough steady forces (F L > Ft), 
a VL's principal vector usually aligns with the direction 
of motion |n|. Though the details of this phenomenon 
are not completely understood, Schmid and Hauger have 
argued that such an orientation is preferred because it 
minimizes power dissipation JTo| | . They focused on the 
vortex motion at high velocities so that the effect of the 
pinning potential can be handled as a small perturbation. 

It is less clear what would happen if the effect of pin- 
ning were more important and the plastic distortions or 
the elastic energy involved during the vortex motion in- 
creased. The perturbative solution would fail to correctly 
describe the vortex behavior and the elastic energy of the 
VL might play a preponderant role. Indeed, in our sim- 
ulations we find that, during the dynamical reordering 
induced by symmetrical forces, a VL's principal vector 
is aligned in most cases perpendicular to the direction 
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FIG. 2: Dynamical ordering, (a) Defect concentration, ridef, 
average absolute- velocity, (V), and average quadratic dis- 
placements of vortices per cycle, (A 2 a;jv) and {A 2 vn} vs. the 
number of cycles, N. (b) n^e/ as a function of F L P/4 after 
100 cycles (filled circles: F L fixed, open circles: P fixed) . The 
point F L P/A = 1 (gray) corresponds to F L = O.OII8/0 and 
is common to both curves. The shaded area marks the initial 
disordered state. The lines are a guide to the eye. 



of motion (see Fig. [l]). In order to show that this phe- 
nomenon is not due to the boundary conditions or to a 
particular distribution of the randomly located pinning 
centers, we performed a series of simulations where we 
slowly cooled the VL from T ~ T m Q down to T = 
while the VL was subjected to the square force in differ- 
ent fixed directions. The temperature was decreased in 
small steps after each cycle. After 600 cycles, T = 0. The 
resulting structure factor of the VL at T — is shown in 
Fig. ||a) for three different directions of the applied force 
(indicated by the white line) . The distribution of defects 
is the same for the three cases. The clear tendency of the 
lattice to choose a transverse orientation (note that the 
structure factor is rotated 90° in relation to the real lat- 
tice) strongly suggests that this phenomenon is not the 
result of either the boundary conditions or an orientation 
favored by the particular distribution of defects. In fact, 
the preferred orientation of the VL can be understood 
with a simple energetic argument [S 14|. In our calcu- 
lations, the vortex motion tends to be disordered and 



neighboring vortices move at different average velocities. 
Some of the vortices can be pinned and might not even 
move during a complete oscillation. The oscillatory force 
produces a back-and-forth movement of vortices and a 
consequent repeated interaction between neighbors. If 
the overlap of vortices is avoided, the average elastic en- 
ergy during one period of the oscillation will be small. 
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FIG. 3: (a) Structure factors of the vortex configuration after 
"annealing" the sample from T ~ T m to T = with a square 
force applied in the directions indicated by the correspond- 
ing white line at the right-bottom corner of each image. It is 
seen that the VL orients with one of its principal vectors per- 
pendicular to the direction of the force, (b) Restitutive force 
considering one vortex moving 5x relative to its neighbors. 
We consider two possible orientations of the VL. 

This favors an ordered structure. Moreover, the most 
favorable condition will be the one that minimizes the 
elastic energy. Because of that, we believe that the elas- 
tic energy is a good quantity to estimate the stability of 
the different orientations of the VL. To this end, we con- 
sider one single vortex moving a distance Sx relative to its 
neighbors, or vice versa, and we calculate the restitutive 
force for different orientations [Fig. 3b)]. From this cal- 
culation, it is straightforwardly seen that the lower elastic 
energy will correspond to the lattice with the principal 
axis transversal to the direction of motion, as observed in 
simulations. The same conclusion is drawn if we consider 
the distortions that are introduced in a nonrigid lattice. 

Similar arguments have been used to explain the col- 
lective motion of colloidal particles in oscillatory shear 
flow Jl5| . The main and very important qualitative dif- 
ference between the colloidal and vortex systems is the 
presence of pinning centers in the latter. However, the 
pinning centers favor an oscillatory shear even when the 
applied force is uniform and, as a consequence, the un- 
derlying mechanism of ordering is probably the same in 
both systems. 

In conclusion, we performed numerical simulations on 
the dynamics of the vortex lattice that emphasize the in- 
trinsic properties of the oscillatory motion. We show that 
the vortex lattice may order when forced to move with 
an oscillating force, even when the motion that produces 
a steady force with the same absolute value is plastic and 
disordered. As the vortex lattice orders, its mobility in- 
creases and it moves in an increasingly coherent way. The 
number of defects in the lattice and its mobility vary as 



the logarithm of the number of cycles of the oscillating 
force. The lattice orders because an ordered structure 
avoids the overlap of neighboring vortices. For the same 
reason, the lattice orients with one of its principal vectors 
perpendicular to the direction of the force. 
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